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I. INTRODUCTION 

The CMBR gives a strong observational basis for the 
assumption of homogeneity and isotropy of our universe. 
Indeed, the deviation of universe from a homogenous 
and isotropic space is ~ 10~ 5 at the stage of decou- 
pling. Thus our universe is well approximated by the 
Friedmann-Lemaitre model before this stage. On the 
other hand the present universe is highly inhomogeneous 
on small scales. The model consisting of the spherical 
shells of dust with the common center, where the rest 
observer is situated, was considered in P, to study 
how the small scale inhomogeneities affect the global dy- 
namics when averaged on larger scales. The gravitational 
field in the region between the shells is described due 
to Birkhoff's theorem by the standard Schwarzschild ge- 
ometry and Israel's Q matching conditions are used to 
derive the equations of motion of the dust shell. The 
dynamics of this model is exactly integrable that allows 
one to make the analytical evaluations. On the other 
hand in the last two decades Milgrom's modified Newto- 
nian dynamics (MOND) paradigm Q has gained recog- 
nition as a successful scheme for unifying much of extra- 
galactic dynamics phenomenology without invoking dark 
matter By presuming the Birkhoff's law the relation 
between modified gravity and corresponding cosmology 
was established in |fg. In this framework the viability 
of MOND cosmology was considered in ■ In this brief 
note we consider the dust shell universe under assump- 
tion of modified Schwarzschild solution. Throughout this 
paper we put G = c = 1. 

Let us begin by considering of a single dust shell under 
assumption of modified Schwarzschild solution given by 



ds 2 = f^dr 2 + r 2 (d0 2 + sin 2 9d(p 2 ) - fdt 2 
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where / = 1 — 2g(m/r) with some function g differing 
from the standard Schwarzschild solution g(x) — x. Us- 
ing Israel's matching conditions [j| the equation of mo- 
tion of a spherical shell of dust takes the form 
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where the prime stands for the derivative with respect to 
the whole argument and over-dot denotes derivative with 
respect to the proper time r. The relation between t and 
r is 



dt 



= [l-2g(m/R)Y 
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The first integral of eq. 10 is given by 



(1 + i? 2 ) 1 / 2 
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where a is the constant of integration characterizing the 
binding energy of the shell. From cq.(0J one sees that 
in the case of standard Schwarzschild solution expanding 
shell with positive binding energy corresponding to a < 1 
falls back after reaching a finite maximal radius. In the 
case of MOND the gravity at large distances behaves as 
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or if accommodated to the cosmological data 

where Ho denotes the present value of the Hubble con- 
stant. So, in the case of MOND the dust shell can tunnel 
through the potential barrier to the regime of unbounded 
expansion. 

Let us now consider the dust shell universe which rep- 
resents the spherically symmetric shells of dust with a 
common center where the rest observer is situated. De- 
noting the mass and radius of the ith shell by m» and Ri 
the equation of motion of the ith. shell takes the form 



R ^Vl + R 2 + ^Jl + R 2 -2g{^j | 
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The first integral of this equation can be written as 



l + Ri = af + g 
1 
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Recalling that in the standard case the matter density p, 
is defined as 0, Q 



mi— i + fTH 



where Xj denotes the initial radius of the zth shell, we 
find it natural to define the matter density in the case of 
modified gravity as follows 
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If one sets the parameters containing in eq.0 as follows 

Pi = p , a? - 1 = -kx\ 
the expansion law of the dust shell can be written as 
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The last radiation like term in eq. Q can be neglected in 
the large N limit 0, Q • Namely, defining 



X; = iAx , NAx 
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where iV denotes the number of shells one finds 



m 2n - m 2n -i = 



m2n+i ~ m 2r , 



n z 108 
W\8^p ' 
6n 2 + 6n + 1 / 3 
87rp 
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So that for large enough values of N the difference 
g (rrii/Ri) — g (m,_i/-Ri) becomes quite small that makes 
it possible to omit this term in the large TV limit. Ne- 
glecting the last radiation like term and introducing a 
new function A 



J 9 ^ 

from eq.ljo]) one finds 
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which after redefinition Xi/Ri — ► scaZe factor determines 
the cosmology with an arbitrary curvature term corre- 
sponding to the modified gravity 0, ■ 

To summarize we presented exact solution of dust shell 
universe dynamics in the case of modified gravity and by 
means of this solution derived the equation governing 
the cosmological evolution. Using the present solution 
one can continue in the spirit of papers 0, Q to study 
an impact of small-scale inhomogeneities on the modi- 
fied cosmology. However, one has to keep in mind that 
this model is highly idealized. In general local inhomo- 
geneities strongly affect the light propagation, giving rise 
to dispersions in the observed distance-redshift relation. 
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